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ABSTRACT 

Nucleation rates in condensed systems are frequently not at their steady state 
values. Such time dependent (or transient) nucleation is most clearly observed in 
devitrification studies of metallic and silicate glasses. The origin of transient 
nucleation and its role in the formation and stability of desired phases and 
microstructures are discussed. Numerical models of nucleation in isothermal and 
nonisothermal situations, based on the coupled differential equations describing cluster 
evolution within the classical theory, are presented. The importance of transient 
nucleation in glass formation and crystallization is discussed. 

I. INTRODUCTION 

Since most first order phase transformations occur by nucleation and growth, 
the prediction of resulting phases and microstructures requires a precise knowledge of 
the rates of these processes. In some cases, such as for many silicate glasses, it is 
possible to measure the nucleation rate directly by counting the number of nuclei per 
unit volume, N v that appear with annealing. N v is related to the nucleation rate, I, by 

t 

N v = f / dt . (1) 

o 


Assuming that the nucleation rate is time-independent, given by the steady state 
nucleation rate, I s , a plot of N v vs the time, t, should then be fit by a straight line with 
slope equal to P. Nucleation studies in many glasses, however, find a nonlinear 
behavior in this plot, generally indicating an initially low nucleation rate that 
approaches the steady state value only after long annealing times. Such transient 
nucleation is common, having been reported in metallic glasses, glass-ceramic 
materials, enamel, and silicate glasses. It is important in many phase 
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transformations, including nucleation of new crystalline phases within existing 
crystalline phases, liquid phase separation, crystallization of undercooled liquids, 
devitrification of glasses, and glass formation. 

In this article, I will discuss the origins of transient nucleation and focus on 
the results of computer studies that directly incorporate transient nucleation to better 
understand phase transformations. Many phenomenon are incorrectly described if a 
steady state nucleation rate is assumed. We show that by using simple numeric 
models, experimental kinetic data can be analyzed in new ways to extract the time 
dependent nucleation rate, I(t), to verify the kinetic model assumed in the classical 
theory, and to investigate glass formation. The interested reader is referred to 
references 1-3 for more comprehensive reviews. 

EL CLASSICAL THEORY - STEADY STATE NUCLEATION 

For simplicity, spherical clusters with a shaip dividing surface, having the 
same composition as the parent phase, are assumed. Following Gibbs, the reversible 
work of formation of a cluster of size n, constructed firom n monomers, of molecules, 
is: 


A G n = nAG' + (36it) 1/3 v 2/3 n 2 * 3 a , (2) 


where AG' is the Gibbs free energy per molecule of the new phase less that of the 
initial phase, v is the molecular volume and a is the interfacial energy per unit area. 
The surface energy penalty is dominant for small clusters, leading to a maximum in 
AG n at the critical size, n*. 


n* = 32u_oi_ 
3v |A GJ 3 


(3) 


where AG V is the free energy difference per unit volume, AG'/v . Since clusters 
smaller than n* will on average shrink, while those larger than n* will on average 
grow, the nucleation rate is in a crude sense the production of postcritical clusters. 

Clusters are assumed to evolve in size by a series of bimolecular reactions, 
describing cluster growth by growth or shrinkage in one monomer steps. The time 
dependent cluster density, N„, is then determined by solving a system of coupled 
differential equations of the form 


£ ~ ^H-1 V^n,l K + ^n, I + ^n*l, t 1 • (4) 

where k^, and k n are the forward and backward rate constants respectively for a 
cluster of size n. Typically they are calculated assuming that the atomic mobility is 
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identical to that describing bulk diffusion. The nucleation rate at a cluster size n, I^ t , 
is the time-dependent flux of clusters past that size, given by 

7 = N k* - N k ' , (5) 

t nf "n, t K n n-l, t K n*l 

Assuming that the cluster distribution reaches steady state, Becker and Ddring 4 arrived 
at the familiar expression for the steady state nucleation rate: 


I 3 = W exp 


AG„. 

kT 


c 


IAG1 


6nk B Ta* 


f, 

H “ K.(D C(T) 


(6) 


where N 0 is Avagadro’s number, k B is Boltzmann’s constant and N n .(T) is the 
temperature dependent steady state distribution. 

HI. CLASSICAL THEORY - TRANSIENT NUCLEATION 

When a glass is rapidly quenched from the melt, there is insufficient time for 
the cluster distribution to maintain the appropriate steady state distribution. Upon 
reheating, then, the density at the critical size will be artificially low, giving a 
nucleation rate lower than expected. With time, the distribution will evolve to the 
steady state distribution appropriate to the annealing temperature, causing the 
nucleation rate to rise. As will be seen in subsequent sections, this process can lead to 
complex nucleation behavior with time. 

Assuming that n is continuous, the coupled set of linear differential equations 
can be written as a single partial differential equation. 


dN n,' _ a 

K K — 

Ml 

dt dn 

" B dn 

i * : 


where N* n is the equilibrium density of clusters of size n. Most expressions for I(t) are 
obtained by solving Eq. (7), making different approximations for the rate constants and 
AG n (see ref. 1). All assume that the relaxation to I s can be characterized by a single 
time constant, the transient time. Experimental data and a computer calculation by 
Kelton et al. 5 (discussion following) suggest that the best of these solutions is the one 
obtained by Kashchiev, 6 
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( 8 ) 


/ = l 3 

«v 


1 + 2 £ (-1)" exp|-^ 

^ F" 


» = 1 


where I s is the steady state nucleation rate, and t k is the transient time, 

-24* Jn* 

x K = ^ . 

rtV.AG' 


(9) 


To compare this with experimental data, however, a proper account must be made of 
the additional time required for a cluster to grow from n* at the nucleating temperature 
to n* at the growth temperature. 1 

Kelton et al solved Eq. (4) by a numerical method to obtain the transient 
nucleation behavior, avoiding the usual restrictive assumptions. 5 Dividing the time 
into a large number of intervals, St, the number of clusters of size n at the end of each 
interval is calculated by 


N. 


V + 5/ 


= N. 




dN 

6 1 — ^ 
dt 


( 10 ) 


Values for I(t) were computed directly from Eq. (5). 

Assuming parameters appropriate to lithium disilicate, good agreement was 
found between the computer generated values for I(t) at n* and those predicted by Eq. 
(8). Recent theoretical treatments by Wu 7 and Shneidman et al„* however, suggest 
that this agreement is fortuitous because of the particular parameters chosen. 
Shneidman argues that the cluster distribution evolution cannot be described by a 
single relaxation time, making most existing analytical treatments invalid. This is 
supported by recent studies by Kelton and Greer, of the effects of preannealing 
treatments on nucleation rates ( cf Sec. V.A). 

Due to its simplicity and more accurate solution of the rate equations, the 
numerical calculation is easily extended to include nucleation under nonisothermal 
conditions, 9 heterogeneous nucleation, 10 and compositional effects. 11 

IV. NONISOTHERMAL NUCLEATION AND CRYSTALLIZATION 

It is important to assess the importance of transient nucleation in 
nonisothermal phase transformations, such as are studied by differential scanning 
calorimetry (DSC) or differential thermal analysis (DTA) for example, and for glass 
formation. To extend the rate kinetics to large clusters, it is assumed that for 
sufficiently large n, N^,= N n+U , giving an average growth rate for a cluster of radius r 
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( 11 ) 



where C is a constant that is generally between 1 and 10, depending on the density of 
the cluster distribution at the point where the extrapolation is begun. Calculated 
growth velocities are in excellent agreement with experiment. 12 

To model a nonisothermal phase transformation, the continuous scan is divided 
into a series of isothermal anneals of size 8T and duration 5t = 8T/S, where S is the 
scan rate. At the end of each interval, the total number of nuclei generated during the 
interval is calculated and stored, and the sizes of the nuclei generated in previous 
anneals are calculated using the average growth rate given in Eq. (11). Assuming no 
impingement, the total volume is computed as 


V, = 



( 12 ) 


where r u is the time dependent radius of each nuclei formed in the previous 
temperature intervals and N; is the corresponding density. Taking V 0 as the sample 
volume, the actual volume transformed is obtained by the Avrami analysis 


a^r) = 1 - exp 



(13) 


DSC and DTA experiments measure the rate of energy production (dH/dt), 
which is proportional to the rate of volume fraction transformed. This is easily 
calculated as a function of temperature as 


DSC signal = 


dH 

dt 


x(T i + 67) - x(T ) 
6 1 


(14) 


V. EXAMPLES 

Three examples are provided to illustrate the utility of these numerical 
approaches. Lithium disilicate is used as a model glass in most cases since it is known 
to nucleate crystals of the same composition by a homogeneous mechanism. The 
parameters determined by Kelton and Greer from quantitative fits to the steady state 
and transient data are assumed. 13 
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A. VALIDITY OF THE KINETIC MODEL IN CLASSICAL NUCLEATION 

THEORY 

Definitive tests of the classical theory of nucleation do not exist, partially due 
to the uncertainty in the values for the interfacial energy, a, and the atomic mobility. 
Most experimental and theoretical attention has focused on the thermodynamic aspects 
of the theory. The kinetic model, common to many theories, however, has received 
little attention. 

Kalinina et al. u first noticed that a proper annealing schedule for lithium 
disilicate glass causes the nucleation rate to go through a maximum with time, 
corresponding to a complicated response of the cluster population to the annealing 
history. Within the classical theory, this behavior should be predictable, providing a 
critical test of the cluster dynamics assumed in the kinetic model for nucleation theory. 
Using the numerical model discussed previously, a and the atomic mobility were 
determined by consistently matching the measured temperature-dependent values for F 
and x. Holding all parameters fixed, then, the number of nuclei, N v , produced as a 
function of annealing were computed and compared with experimental data. 13 

Figure 1 shows N„ as a function of annealing time, measured by Kalinina et al. 
by annealing at 758K and growing at 899K, following preanneals at 713K, 724K, and 
756K. 



Figure 1 The number of nuclei in lithium disilicate glass as a function of annealing time at 
758K after annealing to steady state at 713K, 724K and 746K. (from ref. 13). 

Unlike the case for an as-quenched glass the nucleation rate is too high initially, 
approaching the lower steady state value with time. This demonstrates that the initial 
concentration of clusters of size n* was too large, reflecting the lower temperature 
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distribution induced by the preannealing treatment. The solid lines are the predictions 
from our numerical treatment, showing excellent agreement with the 724K and 75 6K 
data. The apparent disagreement with 713K is likely due to an error in the steady 
state value originally reported at that temperature; the dashed line is the prediction 
based on a value reported later. In agreement with experimental data, calculations 
for lower temperature preanneals predicted a peak in the nucleation rate with annealing 
temperature. Investigations demonstrated that this corresponds to a complicated 
evolution of the cluster distribution that could not be described by a single relaxation 
time, in agreement with the theoretical arguments of Wu 7 and Shneidman. 8 Our 
computer generated results demonstrate for the first time that nucleation proceeds as if 
there were an actual cluster distribution that evolves with the dynamical model 
assumed in the classical theory. 

B. GLASS FORMATION 

Glass formation requires that the liquid be quenched so fast that there is 
insufficient time for significant nucleation and growth of crystalline phases; a glass is 
typically assumed to have a crystalline volume fraction, x c less than Id 6 . 15 
Calculations of the critical quenching rate for glass formation, Q*, are usually made 
assuming steady state nucleation rates and macroscopic growth velocities. 16 During a 
rapid quench, however, it is likely that the steady state distribution will not be 
maintained, resulting in a nucleation rate that can be significantly lower than expected, 
making glass formation easier and increasing the stability of the glass to 
crystallization. 

By directly simulating the evolution of the cluster distribution under the 
nonisothermal conditions of a quench, Kelton and Greer 9 demonstrated that the 
nucleation rate decreases markedly with increasing quench rate (Fig. 2 left). The 
effect is even more dramatic when parameters are chosen that are appropriate for 
Au 81 Si l9 , a typical metallic glass forming alloy (Fig. 2 right). 



Figure 2 Nucletation rate as a function of quench rate for lithium disilicate (left) and Au^Sq, 
(right) (from ref. 9). 
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The lower nucleation rate results in a striking decrease in the number of nuclei 
generated (Fig. 3.a) and in the volume fraction transformed (Fig. 3.b, calculated 
following Sec. IV) over expectations from steady state calculations in Au 81 Si 19 . Steady 
state calculations predict a value for Q* between 10 7 to 10 8 K/s ( for x c < 10" 4 ); 
transient effects decrease Q* to = 10 s K/s, in agreement with experiment. It is 
therefore likely that only because of transient nucleation is glass formation possible in 
some metallic alloys. In lithium disilicate, the predicted values for Q* are identical, 
whether transient effects are included or not, resulting from the decreasing 
crystallization rate for temperatures below that where the nucleation departs 
significantly from its steady state value. 



Figure 3 Number of nuclei and volume fraction transformed as a function of quench rate for 
steady state and transient nucleation in Aug,Si 19 (from ref. 10). 

C. DSC AND DTA MEASUREMENTS OF NUCLEATION RATES 

Based on an experimental study of lithium disilicate glass, where F is well 
known, Ray and Day 17 suggested that nonisothermal DSC or DTA experiments could 
provide a rapid method for obtaining nucleation data, arguing that the temperature of 
the maximum nucleation rate could be determined by first heating the glass 
isothermally in the DSC or DTA at a sufficiently low temperature to promote crystal 
nucleation and subsequently measuring the temperature of the maximum rate of 
transformation, T p and the peak height at T p , 5T p , upon nonisothermal heating. The 
minimum in a plot of T p (or the maximum in 5Tp) vs the annealing temperature, 
should correspond to the maximum in the nucleation rate. Using the numerical 
methods described in Sec. m and IV the experimental procedures were simulated 
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exactly for lithium disilicate glass. 12 Figure 4. a shows the predicted DSC peak heights 
(for a scan rate of 15°C/min) as a function of the preannealing temperature (assuming 
3 hour preanneals). Values for I ! (T), calculated using the same thermodynamic and 
kinetic parameters are also shown. As observed experimentally by Ray and Day, a 
peak is predicted near the peak in the steady state rate. These two peaks are shifted 
slightly, however. Also in agreement, the peak temperatures and peak widths showed 
a minimum near the peak in the nucleation rate. More recent calculations on 
Na20.2Ca0.3Si0 2 , a glass with a much larger nucleation rate, produce results that are 
similar to those for lithium disilicate, but the peak in 8T p is much weaker (Fig. 4.b). 
Again, the peak in the plot of 5T p vs. annealing time is displaced slightly from the 
peak in the steady state nucleation rate, although as shown, the location of the peak is 
sensitive to the duration of the preanneal. 
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Figure 4 DSC peak heights as a function of preannealing temperature for (a) lithium disilicate 
and (b) Na 2 0.2Ca0.3Si0 2 glasses. Computed steady state nucleation rates are shown for 
comparison. Figure (a) is taken from ref. 12. 


These results qualitatively support the experimental conclusions of Ray and 
Day. They depend critically on transient effects. Predicted results from steady state 
calculations are in poor agreement with experimental data. 


VI. CONCLUSION 


In conclusion, transient nucleation is important for a proper understanding of 
isothermal and nonisothermal phase transformations, it is even essential for glass 
formation in some cases. Nucleation rates are a strong function of annealing history, 
often showing complex behavior that verifies the kinetic model of cluster evolution 
assumed by the classical theory. Using the principles developed in this paper for 
modeling this cluster evolution on readily available microcomputers, it is now possible 
to analyze nonisothermal DSC or DTA measurements more accurately than has been 
possible and to devise new methods for extracting fundamental kinetic parameters such 
as the nucleation rate. 
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